Multiple temperature scales of the periodic Anderson model: 
the slave bosons approach 
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The thermodynamic and transport properties of intermetallic compounds with Ce, Eu, and Yb 
ions are discussed using the periodic Anderson model with an infinite correlation between / electrons. 
At high temperatures, these systems exhibit typical features that can be understood in terms of a 
single impurity Anderson or Kondo model with Kondo scale TV- At low temperatures, the normal 
state is governed by the Fermi liquid (FL) laws with characteristic energy scale To. The slave boson 
solution of the periodic model shows that To and Tk depend not only on the degeneracy and the 
splitting of the / states, the number of c and / electrons, and their coupling, but also on the shape 
of the conduction electrons density of states (c DOS) in the vicinity of the chemical potential. 
We show that the details of the band structure determine the ratio To /Tk and that the crossover 
f"*^ ' between the high- and low-temperature regimes in ordered compounds is system-dependent. A 

sharp peak in the c DOS yields To <T]f and explains the 'slow crossover' observed in YbAl3 or 
YbMgCu4. A minimum in the c DOS yields To >Tk, which leads to the abrupt transition between 
the high- and low-temperature regimes in YbInCu4. In the case of CeCu2Ge2 and CeCu2Si2, where 
To — Tk, the slave boson solution explains the pressure experiments which reveal sharp peaks in 
I . the T 2 coefficient of the electrical resistance, A = p(T)/T 2 , and the residual resistance. These peaks 

5— i ■ are due to the change in the degeneracy of the / states induced by the applied pressure. 

' The FL laws explain also the correlation between the specific heat coefficient 7 = Cv /T and the 

slope of the thermopower a(T)/T, or between 7 and the A coefficient of the resistivity. For JV-fold 
degenerate model, the FL laws explain the deviations from universal value of the Kadowaki- Woods 
ratio, Rkw = A/"/ 2 , and the ratio q — UmiT-^oya/^fT. The renormalization of transport coefficients 
can invalidate the Wiedemann-Franz law and lead to an enhancement of the thermoelectric figure-of- 
merit. We show that the low-temperature response of the periodic Anderson model can be enhanced 
(or reduced) with respect to the predictions based on the single-impurity models that give the same 
O high-temperature behavior. 

o 

I. INTRODUCTION 

Intermetallic compounds containing Cerium, Ytterbium and Europium ions exhibit a number of remarkable phe- 
nomena, like the heavy fermion mass enhancement, valence fluctuations, huge thermopower, spin-charge separation, 
unconventional magnetism and superconductivity, etc. These systems have been studied for several decades and are 
still attracting considerable attention. The initial focus was on dilute alloys with magnetic 3c? and 4/ impurities, 
where the low- and high-temperature behavior is qualitatively different. At high temperatures, the experimental 
data indicate that the / states are localized and weakly coupled to conduction states. The susceptibility is of the 
Curie- Weiss form and entropy is large, as expected of nearly free / ions. At the same time, the logarithmic resistivity 
and large thermopower are typical of the conduction (c) electrons weakly perturbed by the local moments. Such 
a behavior is explained by the perturbative solution of the Anderson or Kondo models, which yield the low-energy 
correlation functions as universal functions of reduced temperature T/Tk, where Tk is the Kondo temperature. At 
low temperatures, the experimental data show that the susceptibility is Pauli like, the specific heat and entropy are 
linear in temperature, and transport coefficients are given by simple powers of reduced temperature T/TV/. Close to 
the ground state, the c and / states seem to form a non-magnetic singlet with the fermionic excitation spectrum and 
energy scale Tk- The overall behavior of dilute alloys with paramagnetic impurities is nearly the same when plotted 
on the T /Tk scale, even though the values of Tk vary by orders of magnitude. Unlike the high-temperature data, 
neither the crossover nor the low-temperature behavior can be understood in terms of the perturbation theory which 
treats the conduction electrons and the local moments as separate entities. The fact that the low-temperature scale 
coincides with the high-temperature one and that the crossover from the weak to the strong coupling regime takes 
place around Tk are the most prominent features of dilute Kondo systems. 

The dilute alloy problem has been explained after several decades of intensive work by the exact solutions of the 
Kondo and the Anderson models. Early calculations considered the models in which the conduction electrons with a 
constant density of states (c DOS) are exchange scattered on a spin-1/2 local moment. The solution was obtained by 
the variety of methods, like the perturbative scaling^, the numerical renormalization group^, the Fermi liquid theorjs^, 
and Bethe Ansatz 5 . These results reveal that the effective coupling between the free fermions and degenerate local 
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moments increases at low temperatures and diverges in the ground state, which explains the breakdown of the 
perturbation theory. However, the simple models cannot provide a quantitative description of dilute Kondo alloys; 
a realistic modeling has to take into account the details of the local states and the band structure of the host, and 
consider additional scattering mechanisms. In Kondo systems with Ce and Yb ions, the splitting of the 4/ states by 
the spin-orbit and the crystal field (CF) effects can change the effective degeneracy of the system, reduce the Kondo 
scale, and lead to the seemingly complicated features, like the multiple peaks in the resistivity or the sign change of the 
thermopower and the Hall effect. A non-constant c DOS leads also to a complicated low-temperature behavior. For 
example, Whithoff and Fradkin 6 have found that the c DOS with a power-law singularity leads to a critical coupling 
J c that separates different ground states. For J > J c , the ground state is of the strong coupling type, while for J < J c 
the renormalized coupling decreases with temperature and the "usual" Kondo screening of local moments does not 
occur—. Nevertheless, despite a conduction electron pseudo-gap, it has been shown that Kondo-lattice screening can 
be more stable than in the single Kondo impurity case^. 

The Kondo problem becomes much more difficult for stoichiometric compounds in which a magnetic ion is present 
in each unit cell. The high-temperature features can still be explained by an effective single-impurity model, if one 
takes into account additional splittings of the / states and/or the fact that c DOS can change rapidly in the vicinity 
of the chemical potential. The perturbation theory provides a consistent picture of the high-temperature data: it 
yields the correct Kondo scale, explains the Curie- Weiss behavior of the susceptibility and the logarithmic decrease 
of the resistivity and thermopower, and accounts for the well resolved CF excitations seen in neutron experiments. 
However, at sufficiently low temperatures the scattering becomes coherent and one finds new features that cannot be 
explained by the single impurity models. The onset of coherence is most clearly seen in the electrical resistivity which 
drops to very small values, and in the optical conductivity which shows the development of a low-frequency Drude 
peak and a small hybridization gap close to the chemical potential. At lowest temperatures, the Fermi liquid laws 
often emerge: the resistivity is quadratic and the thermopower a linear function of temperature; the specific heat and 
the magnetic susceptibility are much enhanced, indicating a large effective electronic mass; the de Haas- van Alphen 
experiments show that / electrons contribute to the Fermi volume. The low-temperature ratios of various correlation 
functions, like the Wilson ratio, x/l or the Kadowaki- Woods ratio, A/7 2 , assume the universal values. Here, 7 and 
X, denote the T — > limit of the specific heat coefficient and the magnetic susceptibility, and A is the coefficient of 
the T 2 term in the resistivity. The universality of these ratios indicates that the ground state properties depend on a 
single energy scale To. However, unlike in dilute alloys, this FL scale To can be much different from Tk ■ 

In this paper we study the periodic systems with 4/ ions and show that the relative magnitude of the Kondo and 
the FL scale depends on the shape of the unperturbed conduction states. The Anderson lattice with an enhanced 
DOS around the chemical potential has To <cTx, which explains the gradual transition between the coherent and 
the incoherent regimes (slow crossover) observed in YbAl3 9 . A pseudo-gap or a reduced c DOS close to the chem- 
ical potential yields Tk "CTo, which explains the abrupt valence-change transition observed in Yb- and Eu-based 
intermetallic compounds, like YblnCu^, EuNi 2 (Sii_ 2: Ge a: )2ii or E^Pdi-zPt^Si^. In the case of CeCu2Ge2 and 
CeCu2Si2, where and T seem to be of the same order of magnitude, we show that the pressure-dependence of 
the A(P) coefficient and the residual resistance is driven by the change in the degeneracy of the / states. 

The paper is organized as follows. Section |TT] explains briefly the slave boson formalism for the Anderson model 
and defines its low- and high-temperature scales To and Tk- This section provides also the relationship between 
the c DOS and the relative magnitude of To and Tk , discusses the effects of the magnetic field, and shows how to 
express the transport coefficients in terms of T/To. In Section IHII we discuss the relevance of these results for the 
experimental data on various intermetallic compounds in which the ratio Tq/Tk can be much smaller or larger than 
one. For the case To ~ Tk, we discuss the pressure anomalies in the A coefficient and the residual resistance induced 
by the changes in the effective degeneracy of the / states. 



The periodic Anderson Model (PAM) Hamiltonian is written in the limit of an infinite correlation between / 
electrons, U = +00, as 



II. SLAVE BOSON SOLUTION OF THE PERIODIC ANDERSON MODEL 



A. The mean-field self consistent equations in the limit U — > 00 
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where Ci a and fi a are annihilation operators for c and / electrons, i is the site index, a =|, j is the spin component, h 
is an external magnetic field, and g c and gf the Lande factors for c and / orbitals, respectively. The c orbital describes 
the conduction band with energy levels ek, where k is the momentum component in the reciprocal space, the localized 
/ orbitals are characterized by energy level Ef and the local hybridization between the two orbitals is specified by 
the matrix element V. The infinite Coulomb repulsion restricts the occupation of the / states to rif < 1 and we use 
the chemical potential p to fix the total electronic occupation per site to n = n c + ra/. The unperturbed c DOS is 
Pq{oS) — l/Jv"X)k S{u — £k), where J\f is the number of lattice sites. In the following, pg will be characterized by a half 
bandwidth D. Typical band-shapes and fillings considered in this work are shown in Fig. [TJ All the energies (except 
the excitation energies) are measured with respect to the center of the c-band (see Fig. [1]). The effective degeneracy 
of the model is determined by the lowest spin-orbit state or the additional crystal field splitting of the / states of the 
Ce and Yb ions. This degeneracy can be changed by temperature, pressure or chemical pressure. Here, we consider 
an effective spin-1/2 model. 




FIG. 1: (Color online) Schematic plot of the non-interacting c DOS po(u>). (a) For a regular c DOS and far from the electronic 
half-filling we find To <Tk. Close to the half-filling, To/Tx depends on the shape of the c DOS around the renormalized 
chemical potential /j,: (b) po(u>) is nearly constant for u) ~ p and To ~Tx. (c) \i is close to a minimum of po(w) and To 3>Tk- 
Here, po f» \i is the chemical potential of n c — n— 1 non- interacting c electrons ("small Fermi surface"), and pl is the chemical 
potential of n c = n non-interacting c electrons ("large Fermi surface"). 



The model defined by Eq.(JTJ) is solved for an arbitrary c DOS by the slave boson mean-field (MF) approximatio: 
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which represents the / states by the product of spinless bosons, &T, and auxiliary fermions of spin a =|, j, dj . By 
definition, these auxiliary particle operators create the local / states with no electron and one er-spin electron, 
respectively; \0)f — ► b\ and \a), — ► d\ a |0), r The electronic operators are related to the auxiliary operators as 
fie = b\di a , which leads to the local identities f\ a jia — d\ a d ia . The doubly occupied state |TI)f is forbidden in the 
limit U — ► co. The anti-commutation relations for / operators as well as the physical Hilbert space are recovered 
by enforcing the local constraints b\bi + ^2 a d\ a di a — 1. The latter are satisfied by introducing the time-dependent 
Lagrange multipliers A;(t). Assuming g c << gf we set g c = and gf = 1, replace hgf — > h, and rewrite the PAM 
Hamiltonian (fT]) in terms of auxiliary fermionic and bosonic operators, 
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Since the slave boson Hamiltonian {2j is invariant under the local gauge transformation bi — > bie l6i , di a — > di a e ' 1 , 
we choose the gauge such that the bosonic fields are real, 6,* = b\ = r,. Finally, we make the MF approximation, 
assuming that the boson fields and the Lagrange multipliers are homogeneous and spin-independent constants, = r 
and Xi = A. Within this MF approach, which is exact in the limit of a large number of spin components, the 
Hamiltonian © becomes quadratic, 



TjMF \ " 
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rV 



Ej-X-n 
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(3) 



where 



Aj 



X + h 
X-h 



(4) 
(5) 



In the presence of the magnetic field h, the Lagrange multipliers are shifted by ±/i but the up- and down-spin states 
remain decoupled. The self-consistent solution is obtained by minimizing the free energy [3!F{r, A) = —lnTre~@ H 
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2r(E f - A - M ) - E< C LA CT + tCcta) , (6) 



with respect to r and A. ^From dJ-(r, X)/dr = and dF(r, X)/dX — we obtain 

V ■ 

= i-^£<4a*>> (7) 

kcr 

where (• ■ •} is the thermal average with respect to the MF Hamiltonian Q- The total electron occupation per site is 



n = Jf ^( d L d i<y + c \a c io) = n f + n c ■ (8) 

Here, n f = l/N"52 i(T (fi a fia) = l/A^X^R 1 ^^) and n = = 1 /^J2ia( c L c ^)^ are tne average occupations per site 
for / and c orbitals, respectively. Using the relationship between the slave-boson amplitude and the total electronic 
occupation, 

n f = 1 - r 2 , (9) 
n c = n — nt ■ (10) 

the self-consistency equations ((6|), ([7|), and ([8|) can be written for r ^ as, 

A + n - E 



V 2 £ = ^E]^ pScH«fH^, (ii) 

/+oo 
p° d {u)n F {w)(ko , (12) 
-oc 

/+oo 
p^wJn^Hdw, (13) 
-CxD 



where uf(u)) = 1/[1 + e 73 "] is the Fermi function and p^ c , p a c the spectral densities of the local single-particle 
Green's functions of a given spin component. The Green's functions are defined in the usual way as thermal averages 
of the (imaginary) time-ordered products of the appropriate creation and annihilation operators. For the quadratic 
MF Hamiltonian © , the Green's functions for r ^ are easily calculated by the equations of motion, which yields^ 

= --lm G2» = -^t-pSM , (14) 

7T U) A/j 



p° d {u) = --lm G%(u) = p° c (w) , (15) 



1 / r 2 V 2 \ 
p °(u) = __i TO G -( w ) = p (u; + p . (16) 

7T \ 0J - A CT ) 

In what follows we analyze the MF slave boson solution and discuss the effects of the band-structure. 



B. The high-temperature solution 

At high enough temperatures, the self consistency equations ([6]) - (fTI)]) . have only a trivial r = A = solution, i.e., 
the effective MF hybridization in Eq. ^ vanishes, rV — 0. This solution describes a system of decoupled / and c 
electrons. For h = and Ef < fi, each lattice site is occupied by a single / electron of spin a and there are n c = n — 1 
conduction electrons (per site) with the c DOS given by po(e). The Fermi surface encloses n c = n — 1 points in the 
k-space, i.e., the r = state occupies a "small" Fermi volume that includes "light" c electrons but not the / states. 
The magnetic susceptibility is Curie-like, provided the direct effect of h on the conduction electrons is neglected. (The 
Pauli-like susceptibility of c electrons is negligible with respect to the Curie contribution of the / states.) As long as 
the c states are highly degenerate, T D, their entropy is much smaller than the entropy of the localized / states, 
and the overall entropy per site is approximately given by S = In 2. In the presence of a large magnetic field, the 
degeneracy of the / states is lifted and the system acquires additional Zeeman energy. 
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Even though the trivial r = A = solution does not provide a quantitative description of the PAM at high 
temperatures, it captures the essential qualitative point: it represents the whole system in terms of two well-defined 
but separated sub-systems. A more realistic approach would take into account the small coupling between the c and / 
electrons and treat it as a perturbation. This would reduce the average / occupation to rif < 1, give the Curie- Weiss 
rather than the Curie susceptibility, and obtain the transport properties from the scattering of "light" c electrons 
on the / ions. The resistivity calculated in such a way has logarithmic corrections to the high-temperature spin 
disorder limit. However, both the trivial slave boson solution and the perturbative one break down at sufficiently low 
temperatures. 

Remarkably, the Kondo scale Tk defined by the high-temperature perturbation expansion agrees with the char- 
acteristic temperature at which the non-trivial solution of he slave boson equations emerges. In what follows we 
analyze the non-trivial solution of the self-consistent equations © - (|10j) for the spin-1/2 model and show that it 
captures the main features of the experimental data on Ce, Eu and Yb intermetallic compounds at low temperatures. 
The generalization to an arbitrary SU(N) symmetric Anderson model or the model with the CF split / states, is 
straightforward (see Appendix). 



C. Kondo temperature Tk 



The non-trivial r =/= solution of the slave boson equations is found below some critical temperature which defines 
the Kondo scale Tk of the periodic Anderson model. For a given set of parameters and total occupation n, the Kondo 
scale is obtained from the r, A — > limit of the equations ([6]) — ([8j). This gives, for h = 0, 

Po(u) + p) tanh (ui /2T K ) ^ ^ 



J(fJ-) J-oo U 

/+oo 
dw 8{uj - A) n F {uS) = 1 , (18) 
-oo 

/+oo 
dcu Po(uj + p) np(uj) , (19) 
-oo 

where we introduced the Kondo coupling constant, 

2V 2 

J (p) = -. (20) 

H-Ef 

Assuming that Tk vanishes continuously and taking the limit Tk — > yields the critical coupling 

2 f + °° , p (u> + p) 



J, (21) 

For a regular DOS, po(p) ^ 0, the r.h.s. of Eq. ([2Tj) diverges logarithmically, such that J c = and any finite coupling 
leads to Tk > 0. At Tk, the high-temperature r = phase with large paramagnetic entropy is destabilized by a 
transition to the low-entropy Fermi liquid phase. The Kondo scale given by Eq. (|17|) coincides with the solution of 
the scaling equations for the SU (2) single-impurity Kondo model^. 

The physical interpretation of the emergence of the non-trivial slave boson solution is made by the analogy with 
the single-impurity case. We assume that the high-temperature solution of the lattice describes a system of localized 
/ and c electrons with small FS, and that the correction to the r — solution can be found by the perturbation 
theory in terms of J{p). For T > Tk the hybridization energy is small with respect to the entropic contribution to 
the free energy due to the degenerate / states. Thus, the total free energy of the system, T ' = E — TS, is minimized 
at high temperatures by the paramagnetic configuration in which the / states are very weakly coupled to the c states. 
For T < Tk, the entropic contribution is reduced below the hybridization energy and the paramagnetic configuration 
becomes thermodynamically unstable. At Tk, the crossover between the local moment and a local singlet state takes 
place. 

In the case of an unperturbed c DOS with a pseudo-gap at p,, po(p^) = 0, equations (TTj) - (fT5|) yield the solution with 
finite J c which separates two different low-temperature regimes. For J > J c the r ^ solution emerges at temperature 
Tk which is the same as in the po(p) ^ case. But for J < J c the r = solution persists down to T = 0, such that the 
paramagnetic entropy is not removed by Kondo scattering. If the coupling is tuned by pressure or doping, a quantum 
phase transition can be induced at the critical value J c (or V c ). In the case of c DOS with a gap A g <C D around 
the Fermi level we find J c oc D/ln[D / A g ]. A pseudo-gap centered at the chemical potential p and characterized by a 
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single energy scale D gives J c oc D. The power-law singularity po(w) — Ro\w\ n gives J c = 2i]D/(r] + 1), which is the 
Whithoff and Fradkin result^. (The constant Rq follows from the normalization condition, J D D Pq(lu) = 1.) 



Se(T) 




FIG. 2: (Color online) Schematic plot of the electronic contribution to the entropy S e - as a function of the normalized 
temperature T/Tk for three cases: To >Tx (red dash doted line), To ~Tk (blue solid line), and To << Tk (green dashed 
line). The doted lines indicate the linear Fermi liquid regime S e -(T) — T/Tq, with a rescaled slope Tk/To. For T>Tk the 
three curves are identical, reflecting the linear contribution from the conduction band S e - (T) — In 2 + T/D. The black dot 
refers to a standard determination of Tk from the electronic entropy: S e -(Tk) = a; In 2. On this schematic plot, we used 
x — 1 , which gives Tk as defined by the slave boson (r = 0) solution. Experimentally, where Tk is defined by the crossover, 
one usually takes x = 1/2. 

To find Tk at finite po{p) we use the Sommerfeld expansion of equation (jTTJ) which gives in the Tk « D limit 
the result 

T K = a {D 2 - pl) l l 2 F K e- l ' J ^ p0 ^ , (22) 

where 



Fk — exp 



-doj 



-(D+Ho) 2|w|p (^o) 

u is measured with respect to p, ao = 1.13 is a numerical constant, and po is defined by the integral 

n c n — 1 



(23) 



2 



Mo 

p {uj)du . (24) 

-D 



By definition, po is the Fermi level of n c non- interacting electrons which have a "small" FS. The result given by 
Eq. (|22[) is derived in Rcf. 16 for a Kondo lattice with an analytic DOS. It also holds for po(oj) with an algebraic 
singularity close to po. 

The Kondo temperature Tk, defined in Eq. (fl~7|) characterizes a second order phase transition. We are aware that 
this transition is an artefact of the slave boson MF approximation and that an exact theory would lead to a crossover 
instead. Experimentally, the Kondo temperature Tk is sometime estimated from the resistivity measurements, which 
show a maximum around Tk, or from the specific heat measurements, where Tk is identified as the temperature at 
which the magnetic entropy (per impurity) becomes a substantial fraction of In 2, say S ~ 0.5 In 2 (see FIG. [2]). Note 
that with this definition of Tk from the magnetic entropy we implicitly neglect the collective freezing of the entropy 
which is related to the Ruderman-Kittel-Kasuya-Yosida interaction. 

The fact that the Anderson lattice and the single impurity Anderson model have the same Kondo scale indicates that 
in both cases Tk separates the paramagnetic high-entropy phase from the low-entropy phase in which the conduction 
electrons start forming an incoherent screening clouds which reduce the local / moment in each unit cell (see FIG.[lJ. 
At temperatures at which the c and / states form a coherent band, and local screening clouds become correlated, the 
hybridization V cannot be considered as a perturbation. For T <C Tk, the Hamiltonian has to be diagonalized by 
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non-perturbative methods and the slave boson solution provides a reasonable description of the renormalized ground 
state. Close to the ground state, the low-energy excitations of a periodic system are Bloch waves and we expect them 
to be characterized by a FL scale T . The question is, how is T related to T#. 



D. The Fermi liquid scale To 



The emergence of a strongly coupled / — c fluid at does not imply that Kondo scale characterizes the behavior 
of such a fluid close to the ground state. Electrons described by Eq. ([3]) form at T — a coherent Fermi liquid with 
an energy scale To, which is defined in the absence of a magnetic field as, 



To 



1 



1 



pHp) pKp) 



(25) 



XI 



are given by simple powers of reduced temperature T/T . The slave boson result for T is computed from Eqs. 
IT6| at T = 0, which gives 



where p a {p) = p%{p) + P'fip) 1S the renormalized density of electronic states at the Fermi level. The scale T 
is relevant for the T — > properties of the periodic Anderson model; it determines the static spin susceptibility, 
(T = 0) ~ 1/Tq, the specific heat coefficient 7 = Cv/T ~ l/T), and appears in the transport coefficients^ which 

CD 

(26) 
(27) 



Ef-X-fi 
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n c + rif 



u> — A 



Co w + f 



r 2 V 2 
uj — A 



dui , 



2ta2 



r*V 



(lj - A) 2 



PO [ LU + [I - - - - I (U' = 



po{uj)du} 



where po{oj) 
have 



is the unperturbed c DOS and ph the chemical potential of n = n c 



p, L = fi + 



r 2y2 



A 



non-interacting electrons. We 



(28) 



which defines the shift in the chemical potential due to the enlargement of the Fermi volume of hybridized (r 5^ 0) 
system with n particles with respect to the Fermi volume of the non-interacting (r = 0) band with n — 1 particles, 
A/i = fiL — M = r 2 V 2 /A . (This interpretation of Ap neglects the width of the free-electron distribution function 
at temperature with respect to D and assumes p ~ po, which holds for n/ ~ 1.) The low-temperature FS is 
"enlarged" with respect to the high-temperature one, because it accommodates nj additional / electrons. Equations 
(113, C5||, and (EHD yield 



rV J 



Po(pl) 



(29) 



which shows that p(p) and To depend on the shape of the unperturbed c DOS and the total number of particles. 
In the FL regime, the entropy behaves as T/Tq, the susceptibility is constant and transport coefficients are given by 
simple power laws of T/Tq. The schematic phase diagram of a system with Tk 7^ Tq is shown schematically in Fig. [3] 

The effect of the c DOS can be computed analytically in the V « D limit, assuming that p is close to pq. Since pl 
corresponds to n non- interacting electrons ("large" FS) and po to n c ~ n — 1 electrons ("small" FS), we approximate 
Ap cx D (see FIG.Q}. In the limit rV/D <C 1, the MF Eqs. (EBJ) and give 



E f - A 



2V 2 
n c + rif 



Po(p) In 



r 2 V 2 



(D + p)Ap 
p (p + uj)cLj 



Am Po(p 



Po(p) 



diJ 



where we neglected the corrections of order (rV/D) 2 . Solving Eqs. (|3"01 |3T|) for A yields 

r 2y2 



A 



Ap 



{D + p ) F e -V^o)po(Mo) 



(30) 
(31) 

(32) 
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FIG. 3: (Color online) Schematic phase diagram of the PAM. The reduced temperature T/Tk is plotted versus To/TV, which 
is considered as a tunable parameter that can vary with the electronic filling, magnetic field and/or the shape of the non 
interacting DOS. The crossover regime can be non universal or non Fermi liquid. 



where 



Fn 



exp 



-{D+no) 



Mpo(mo) 



du> 



(33) 



The scale To of a system of particles described by the periodic Anderson model close to the ground state is obtained 
from equations (|25I29I32[) as, 



To 



r-21/2 



1 



D + PO F 

(A^) 2 p {p L ) A/i Po(pl) 



-1/J(mo)po(mo) 



(34) 



The scales To and Tk have the same exponential dependence on the coupling constant J(po) but their pre-factors 
are not affected by po(u)) in the same way and can differ considerably. The ratio T /Tk which is constant for a given 
set of parameters can be changed by applying pressure or magnetic field. 

At temperatures low with respect to Tfl — mm {To, Tk} the system behaves as a Fermi liquid (provided we are not 
too close to the Kondo insulating state) but the crossover from the high-temperature to the low-temperature regime 
proceeds differently for To <^Tk than for To ^T^. This is shown schematically in Figs.[2]and [H If To <CTx, we have 
Tfl = To and, for Tq < T < Tk, the system exhibits an extended non- universal behavior. If the electronic occupation 
n c is neither too small nor too close to half-filling, the MF solution of the slave boson equations^ predicts the entropy 
with a plateau at about (1 — n c ) In 2 (see FIG. [2|), which characterizes 1 — n c unscreened magnetic ions. If To ^> Tk, 
the high-temperature perturbative regime persists all the way down to Tk, where the properties change abruptly and 
the system enters the FL phase characterized by To. Only for To ~ Tk, the lattice system is characterized by a single 
energy scale, as in the single impurity case. 



E. Comparison of the Fermi liquid and the Kondo scales 



The Kondo coupling J(p) defined by Eq. (|20|) is always smaller than the half bandwidth D of the non-interacting c 
DOS, such that Tk and T arc exponentially small due to the factor exp{ — 1/ J(/i )po(Mo)} • The analytic expressions 
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given by Eqs. 



yield the result 



Po 



T K \D - p 



1/2 



F 



a a p a (p L )i^p F K 



(35) 



which does not depend on the Kondo coupling J(p) but, as discussed previously^, varies with the electronic occupation 
and with the shape of the non-interacting c DOS. Electronic filling effects have been discussed using the MF analysis 
of the Kondo lattice^ and the DMFT solution of the periodic Anderson models. In the limit n c — * (n « 1), the 
first factor D + po vanishes and Tq << Tk (see FIG. Q] (a)). Another physically relevant case is found in the limit 
n c ps 1, which corresponds to a Kondo insulator with the vanishingly small po(/ i i)- Equation (|35p yields To >T^- 
and the break-down of the FL laws is neither due to the proliferation of the quasi-particle excitations (controlled by 
To), nor to the thermal destruction of the Kondo screening (controlled by T^), but is due to the proximity of the 
chemical potential to the hybridization gap in the DOS. 

We consider now in more detail the effects due to the shape- variation of po(u>). In order to separate this effect 
from the electronic occupation effects 19 ' 20 , we assume that n c is close to 1 (but not at half-filling exactly, so that the 
system is metallic). The first two factors on the r.h.s. of Eq. (|35l) are then of the order 1, such that To/Tr- sa Fq/Fk- 
Assuming Ap m D — po we obtain from Eqs. (|23[) and (f3"3"| a simple relation 



exp 



° Mo PoiPo + - p (po) 



(D+no) 



2\u\p (po) 



(36) 



which describes the dependence of Tq/Tk on the specific form of po. A constant po gives Tq ~ Tk, which explains the 
T/Tk scaling observed in some heavy fermion compounds. If po is close to a local maximum of po(u>) the integrand 
in Eq. (|36|) is negative in the main part of the integration range, such that Tq < Tk, as found in the systems with 
the 'protracted screening— A sharp spike in Pq(uj) close to po would exponentially reduce T with respect to Tk- 
On the other hand, if ^ is close to a local minimum [see FIG. [1] (b)] one finds T 3> Tk, which could be understood 
by the following intuitive argument. The incoherent Kondo cloud which begins to form at T w Tk involves a few 
conduction states around the Fermi energy po- These states are part of the conduction band with a "small" FS 
(the FS of n c — n — 1 non-interacting electrons). When temperature decreases much below Tk, the local / orbitals 
hybridize with the conduction electrons to form a coherent Fermi liquid which is characterized by a "large" Fermi 
surface (the FS of n non- interacting electrons). Thus, To is affected by all the states between the "small" and the 
"large" FS, as well as some additional holes inside the "small" FSi£. For n c ~ 1 and po close to the minimum of 
po(u>), the low-temperature increase of the Fermi volume leads to the DOS which is much larger than the one used to 
evaluate Tk (see FIG. Q](c)). In that case, the formation of the Fermi liquid is "self-amplified", yielding To ^> Tk- 
We recall that for To ^ Tk the FL regime sets in at temperatures that are the smaller than either To or Tk- 

We are aware that corrections to MF slave boson analysis might occur from a more accurate treatment of the model; 
nevertheless, since the SB approximation is known to be correct at low energy, we expect such corrections to provide 
a similar integral expression, where the l/|w| divergency would be smoothed out at hight frequencies. The analytical 
expression (I36[) still provides a good quantitative estimation of the band shape effect in the vicinity of the chemical 
potential. The relative magnitude of To and T^ is related to the functional form of po(to) around uj = p. 



F. Effect of a magnetic field 



We next consider the slave boson solution in the presence of an external magnetic field which couples to the / orbitals. 
The direct effect on the c electrons is neglected, even though they can be polarized due to the interaction with / 
electrons. In the linear response regime, the local magnetization of the / orbitals, m z {h) = Si ((/-it /it) ~ (/i\/a))> 
is proportional to the applied field h and, in the Kubo formalism^!, the proportionality factor is equal to the local 
static susceptibility, computed in the absence of the field. At T = 0, the magnetization is m z (h) oc h/To, where 
To is defined by Eq. (|25|) Going beyond the linear response, the critical magnetic field h c (T) is defined at a given 
temperature by the transition between the r ^ and r — state. Solving the MF equations (TTTj) - (TT3"]) in the r — » 
limit and using the definition of the Kondo coupling J(p) in Eq. (|20|) , we find for rif = 1 



1 



U3* 



hi 



Po{oj + p) tanh 



2T 



du) 



(37) 



which yields h c ^ for any regular po(oj). Unlike the zero-field result, Eq. (|2"Tj) , which is logarithmically singular and 
gives J c = for any finite po(p), equation (f3"T)) yields J c ^ for any finite field. For h > h c (T), there is only the 
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FIG. 4: (Color online) Schematic phase diagram of the PAM as a function of a magnetic field h. The red solid line indicates 
h c (T) which separates the trivial r — solution from the non-trivial one, r/0. 

trivial r = solution which describes ferromagnetically polarized /-moments decoupled from the conduction band 
(see FIG. HJ. For h < h c (T), the equations have a non-trivial r ^ solution which describes a partially polarized 
Fermi liquid. At T = 0, the weak coupling limit (J <C D) yields the universal relation 



Since «o = 1.13, we have hP c ~ Tk- At finite temperatures, the solution of Eq. (|5T)) with constant c DOS yields the 
critical line 



which separates the trivial solution (the decoupled phase) from the non-trivial one (the FL phase) and holds for any 
J(/i) and filling n. The critical line h c (T) obtained for a constant c DOS is represented schematically in Fig. |U The 
same relation is also found at half-filling, for any c DOS. In general, we expect a 'nearly' universal phase boundary, with 
some small deviations due to the structure of the c DOS and the particle-hole asymmetry. Studying the development 
of the system as a function of temperature at constant field gives the critical temperature Tk{Ji). Changing the field 
at constant temperature gives h c {T). 

The magnetization m z (h) obtained from the slave boson solution at T = is plotted in Fig. [5] as a function of 
reduced magnetic field Ii/Tk for three typical cases: T <Tjc T « Tk, and T » Tk- A non-constant c DOS 
leads to different types of magnetization curves which resemble the temperature dependence of the entropy depicted 
in Fig. [2] (with the T/T^ axis replaced by h/hP c ). For h < min {To, TV}, the system is in the linear response regime, 
shown in FIG. E]by dotted lines; the slope of m z (h) is T^/T , when plotted versus H/Tk- This regime is analogous 
to the FL regime described from the entropy. In higher fields, the behavior of the slave boson solution depends on the 
ratio Tq/Tk- For To -CT^, the T = magnetization is linear for small fields and then saturates. For To < h < Tk, 
the magnetization might have a plateau which signifies a saturated FL (r ^ 0) with a 'large' FS. At very high fields, 
h > Tk ~ h®, the magnetization becomes a universal function of h/h®. In the opposite case, Tk -C T), the low 
field limit gives m z (h) oc h/To and the initial slope appears very small when plotted versus h/h®. Once the field 
exceeds the critical value, hP c ~ Tk, the local moments are unscreened and m z (h) rises rapidly towards the free-ion 
value. Thus, at about h ~ hP c such systems exhibit a meta-magnetic transition from an unpolarized Fermi liquid to 
the polarized spin-lattice (the r=0 phase). Of course, this simple considerations should be corrected for direct and 
indirect effects due to the conducting sea. 



h° c = h c (T = 0) 



a 1 Tk ■ 



(38) 



[h c {T)/hlf + [T/T K ] 2 = 1 , 



(39) 
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FIG. 5: (Color online) Schematic plot of the magnetization m z as a function of reduced magnetic field h/TK for To 3>Tr- (red 
dash doted line), To ~Tj? (blue solid line), and To <CTk (green dashed line). The doted lines indicate the initial slope in the 
linear response regime, where m z (h) = h/To. The black dot refers to the complete polarization of the local / electrons, with 
m z = 1/2, which occurs a the critical field h° c — Tk- 



G. Transport coefficients in the FL regime 



The slave boson Hamiltonian in Eq. ((3|) provides the approximate FL scale To and the low-temperature thermo- 
dynamics but has no relaxation mechanisms that could lead to stationary heat and charge currents. To calculate the 
transport properties of the SU (N) symmetric Anderson model in the T — > limit we use the Fermi liquid theory 22 
which takes into account the quasiparticle (QP) damping and leads to a finite relaxation time. The QP excitations 
of the full model have the same dispersion as the excitations of the MF slave boson model but are restricted to the 
immediate vicinity of the Fermi level. In the u> — > limit, where the imaginary part of the / electron self energy 
£/(w) can be neglected, the QP and the slave boson dispersion assume the same form. The correspondence is ob- 
tained by identifying rV with \fZ] and A with ujf, where ZJ 1 = [1 — d£//dtt>| u= o] is the renormalization factor, 
Qf = [Ef + ReEy(O) — fi]Zf the renormalized position of the / level, and the excitation energies u> are measured in 
both cases with respect to the renormalized chemical potential /i. Unlike the infinitely long-lived MF excitations, 
which are formally defined for u> < D — /i, the QP excitations are defined for Im S/(w) <C w 2 - 

We calculate the transport coefficients of the periodic Anderson model with constant hybridization using the 
fact that the charge and energy current density operators satisfy the Jonson-Mahan theorem 2 ^. This allows us to 
express the charge conductivity by cr(T) — e 2 NL\\, the thermopower by a(T)\e\T = — L12/L11, and the electronic 



contribution to the thermal conductivity by k(T)T 
introduced the transport integrals: 



N(L 2 2 — L\ 2 /Lii). In each of these expressions we have 



L, 



du> 



driF 



m+n — 2 



A(w,T). 



(40) 



where uf(u>) = 1/[1 + cxp(/3w)] is the Fermi-Dirac distribution function and A(w,T) is defined by the Kubo linear 
response theory 21 . At low temperature (—driF/duj) approaches delta function and the main contribution to the 
integrals in Eq. (|40[) comes from the low-energy excitations within the Fermi window, |a>| ^ T. In the u>,T — > limit, 
a straightforward calculation yields for the three-dimensional systems^ 



A(u,T) = -v 2 f p c (w)t{u,T), 



(41) 



where we introduced the unrenormalized velocity t>k = — Vk£k and denoted by vp the average of w 2 . over the 
renormalized Fermi surface of hybridized states. The renormalized c DOS is p c {^) and the relaxation time is 
1/t(oj,T) ~ Im S c (cj), where X c (w) is the self energy of the c electrons which must include the quasiparticle 
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damping. The integrals in Eq. (|40|) are evaluated by the Sommerfeld expansion (for details see RefJi) which yields 
the transport coefficients of the periodic Anderson model as simple powers of reduced temperature T/T . The pre- 
factors of various powers are functions of po(fJ-L), n f, and A/i which depend on the local self energy which is 

difficult to obtain for the excitation energies within the Fermi window, |cj| < T. To avoid this problems, we replace 
To and all other renormalized quantities that appear in the FL expressions by the slave boson MF results. 
The FL result for the electrical resistivity can be written as^ 

^ 9 2 (Ap/n f ) 2 9^(Ap/n f ) 2 Tl ( T V _ 2 

R{T) " N{N-l)^vl {lT) ~ N(N~l)e 2 v 2 F T 2 fej ~ AT ' (42) 

where 7 = (it 2 /6)Np(p) = tt 2 /3Tq and the in the second equality R(T) is expressed on the reduced temperature 
scale T/Tk- The coefficient A — R(T)/T 2 depends not only on specific heat coefficient 7 but on the difference in 
the chemical potentials of unhybridized and hybridized Bloch states, the effective degeneracy of the model, and the 
square of the Fermi velocity averaged over the hybridized FS. This result can be used to explain the pressure- and 
magnetic field-dependence of A, which has been studied experimentally in various systems. 
The FL result for the Seebeck coefficient is given for n/ ~ 1 by the expression 

a(T) = T j-r- = T — ■ (43) 

n/|e| n f \e\T T K 

Since the doubly occupied / states are removed from the Hilbert space, the model is highly asymmetric, and the initial 
slope limj^o a {T)/T never vanishes. In bad metals with a low-carrier concentration, a(T) could be very large. 
The thermal conductivity in the FL regime reads, 

k(T) = Ta(T)£ Q (T) (44) 

which yields in the T — > limit the Wiedemann- Franz (WF) law, k(T) oc Ta(T). However, the usual Lorentz 
number, Lq — 7r 2 /3e 2 , is here replaced by the effective one, 



4)(T) = ^ 



32tt 2 / T x 2 



n 2 VTo 



(45) 



The correction given by the square bracket could lead to substantial deviations from the WF law much below To, 
since the factor multiplying the T 2 term is very large. 



III. DISCUSSION OF THE EXPERIMENTAL DATA 



The slave boson solution of the periodic Anderson model is used in this Section to discuss the effects of the band 
structure on the properties of the intermetallic compounds with 4/ ions. First, we consider the case of the c DOS 
with a maximum in the vicinity of the chemical potential, such that the Fermi liquid scale is much smaller than the 
Kondo scale, To CTk. These results provide a qualitative explanation of the experiments on YbAl3 or YbMgALi. 
The opposite case, To >Tjf, occurs when the chemical potential is close to the minimum of the c DOS; these results 
explain the main experimental features of YbInCu4-like compounds. Finally, for To ~ Tk , which seems to characterize 
CeCu2Si2 and CeCu2Ge2, we discuss the rapid variation of the T 2 coefficient of the electrical resistance and of the 
residual resistance with pressure. The anomalies are related to the pressure-induced change of the effective degeneracy 
of the / states. 



1. The To <C Tk case 



Unlike the experiments on dilute alloys, the overall temperature dependence of the experimental data on YbAl3, 
Ybi-j; Lu-e AI3 9 for x < 0.05, YbMgCu4 and YbCdCu4~, cannot be explained by a single energy scale. In 
these compounds, the high-temperature resistivity is a large and slowly varying function of temperatur e 25 ' 26 , the 
thcrmopowe r 25 ' 26 has a broad peak around 350 K, the magnetic susceptibility and the specific heat have a broad 
maximum at T max ~ 125 K, and the Hall coefficient is typical of a metal in which the c electrons scatter on local 
moments — 122. A similar high-temperature behavior is seen in YbXCu4 (X=Cd, Mg, Zn)2£. The inelastic neutron 
scattering data on YbAl 3 show a broad Lorentzian spectrum centered at about 540 me V 29 ' 30 , which can be understood 
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in terms of the Kondo effect with T# > 500 K 27 . Below 100 K, the electrical resistivity of YbAl3 decreases due to the 
formation of a coherent ground state and for T <T ~ 40 K the resistivity is a parabolic and the thermopower a linear 
function of temperature 25 i 26 i 31 ; indicating a Fermi liquid. For T < To, the zero-field anomalies in x(T) and Cy(T) 
yield an enhanced effective mass^Z, of the order of 1 /To and the neutron scattering dat a 29 ' 30 show a narrow peak at 
about 30 meV, which is due to a hybridization gap. However, if we plot the low-tcmperature transport coefficients on 
a reduced temperature scale T /Tk they appear to be strongly enhanced with respect to the predictions of the single 
impurity Anderson model with Kondo scale Tk — 500 K. Neither the A coefficient of the resistivity nor the slope of 
the thermopower can be explained by the single impurity c calculations that would capture the main features of the 
high-temperature data and give Tk — 500 K. Optical conductivity at 7 K shows a narrow Drude-like response that 
is often found in heavy fermion systems 3 - 2 ^ 3 - 3 - and another mid-infrared peak (MIR) that can be associated with the 
hybridization gap. The optical spectra do not change appreciably for 7 K <T< 40 K, as expected of a system with 
the characteristic energy scale To ~ 50 K. However, the Drude peak broadens and the MIR peak vanishes at higher 
temperatures. The de Haas - van Alphen experiments performed up to h ~ 40 T, show that the effective mass is 
reduced along certain directions in k-space by a factor of 2 without a significant alteration of the shape of the Fermi 
surface 3 ^. The fact that the low-temperature susceptibility anomaly is suppressed in the fields of about 40 T indicates 
that the high-field mass renormalization and the zero-field anomalies below To describe different aspects of the same 
physical state. The low-temperature anomalies are easily destroyed by disorder (they vanish in Yh\~ x Lu x AI3 for 
x > 0.05), which is another indication that they are related to the coherent stated. 

The overall shape of the experimental data and different characteristic energy scales found at high and low tem- 
peratures can be explained by the slave boson solution of the periodic Anderson model close to half-filling. As shown 
in Sees. Ill El and III F[ if the chemical potential is close to the maximum of the unperturbed c DOS, the mean field 
equations yield Tq <C Tk — h a c and give rise to a 'slow crossover' from the LM to the FL phase. Taking To ~ 50 K 
and Tk — 500, as suggested by the experiment, we find that the low-temperature susceptibility and the specific heat 
coefficient are FL-like and much enhanced with respect to the expectations based on the behavior in the incoherent 
regime. The FL results for the T 2 term of the resistivity (see Sec. Ill Gj) explain the enhancement that one finds 
below 50 K when p(T) is plotted on the T/Tk scale. The low-temperature thermopower a(T) cx (Tk/To)T/Tk 
is also enhanced by Tk/Tq with respect to the predictions of the single impurity calculation that reproduce the 
thermopower maximum above 300 K. For temperatures between To and Tk, the magnetization and the magnetic 
entropy are reduced with respect to the free-ion value, which is only recovered for T > Tk- The slave boson solution 
gives m z (h) ~ h/T for h <^ T and m z (h) oc const for h ^> T . Since T <C hP c , the slave boson order parameter is 
finite for h > Tq, i.e., for Tq < h < Tk, the system is a polarized heavy FL with a 'large FS'. Such a behavior explains 
de Haas - van Alphen data which show that the FS does not change much up to the fields of about 40 T. 

The model with the chemical potential close to the maximum in the c DOS describes Ybi- X Lu^A^ for x < 0.5. 
For higher concentrations of Lu there are so many additional holes in the conduction band that the chemical potential 
shifts away from the peak in the c DOS. In that case, To and Tk approach each other and the 'slow' crossover does 
not occur. For x > 0.5, one can expect a crossover from a lattice regime with To =/=Tk to a universal dilute regime 
with Tk —Tq, as discussed in Ref 3 -^. 

The slave boson dispersion defined by the Hamiltonian in Eq. Q explains the Drude and the MIR peaks found at 
7 K in the optical conductivity. However, we hesitate to discuss the temperature dependence of the hybridization gap 
using the mean field results. These results are obtained by enforcing the constraint n/ = 1 only on the average, so 
that the auxiliary fermions are mapped on a free electron gas. The proper redistribution of the spectral weight should 
not neglect the QP damping and should use the solution which is valid at all energy scales. 

2. The To » T K case 

Unlike YbALj or Ybi^^, Lu x AI3 for x < 0.5, the transition between the LM and the FL phase in Ybi_ x Lu x AI3 
for x > 0.7, YbTlCu^^l, YbInCu4 and Ybi_ x Y 2 ,InCu4 for x < 0.5i2 takes place at temperature T„ in an abrupt 
way. In YbInCu4, which we take as a typical example, a first order valence change (VF) transition occurs at ambient 
pressure at temperature T„=40 K. Above T v , the Yb ions are in the 3 + configuration and the magnetic response can 
be understood assuming an independent 4/ hole in each unit cell. The ensuing high-temperature effective moment 
is then close to the free ion value, as observed experimentally. The electrical resistance is very large and nearly 
temperature independent. The Kondo scale deduced from these data is so much smaller than T„, that the / and 
c states are effectively decoupled for T>T„. The magnetic fields up to 40 Tesla do not produce any appreciable 
magnetoresistance which is also an indication of a small Kondo coupling. The large values of the electrical resistance 
and small Hall constant cannot be explained in terms of the spin disorder scattering but can be taken as an evidence 
that the c DOS has a pseudogap or a deep minimum in the immediate vicinity of the chemical potential. When the 
compound is cooled down to temperature T„, the lattice expands and the Yb configuration changes from 3 + to 2.9 + . 
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In the VF phase, the susceptibility and the specific heat coefficientare moderately enhanced (7 ~ 50 mJ/mol K), the 
resistivity is quadratic and the thermopower linear in temperature. The ratio a/T~f is typical of a normal FL2£ but 
the Kadowaki- Woods ratio is anomalously low. The characteristic energy scale in the VF phase is much larger than 
T„ ; the data give To rj 500 K-. The optical conductivity 3 ^, the Hall effect^, and the thermoelectric power— indicate 
that a bad metal with only a few states close to the chemical potential is transformed at T„ into a good metal with 
small 7. That is, the VF transition is accompanied by a major reconstruction of the conduction states. 

An application of pressure shifts T v (P) to lower temperature o 41 ' 42 ' 43 ! 44 without changing the properties of the high- 
temperature state. For T > T V (P), the susceptibility and the specific heat of YbInCu4 can be explained by the crystal 
field (CF) theory of independent / states^ split by Acf — 40 K into an excited quartet and two nearly degenerate 
doublets. This CF scheme agrees with the neutron scattering data 4 ^. The entropy of the LM phase, obtained by 
integrating Cy(T) / r 43 i 45 , decreases from S ~ i?ln8 to S ~ i?ln4 as the system is cooled down to T V (P), as expected 
for two CF quartets without any Kondo screening. This indicates, once again, that the Kondo scale of the LM phase 
is much smaller than T V (P). Pressure stabilizes the paramagnetic phase and at a critical pressure of P c = 2.5 GPa 
the LM persists down to T c = 2.5 K, where the magnetic transition removes the paramagnetic entropy of about 
S ~ i?ln4. The NM R 41 ' 44 and neutron scattering data show that the local moments in the magnetically ordered 
(MO) phase are somewhat smaller than in the LM phase just above Tjv- Since the specific heat coefficient and the A 
coefficient of the resistivity are much larger for P > P c than for P < P c , we speculate that the transition in the MO 
phase is accompanied by an increase of Kondo coupling. As regards the doping, replacing Yb with Y or Lu ion a 42 i 47 
produces similar effects as pressure, i.e., doping stabilizes the LM configuration and shifts T c to lower temperatures. 
In Eu-based intermetallics, the valence-change transition exhibits similar feature a n i 12 as in YbInCu4, except T„ is 
higher (T v > 100 K) and the valence state of Eu ions changes almost completely [from 2 + (f 7 ) to 3 + (/ 6 )]. 

We close this experimental summary by mentioning that the VF transition shifts in an magnetic field to lower 
tempeatures. At ambient pressure, the critical field of h° c ~ 40 Tesla suppresses the VF transition and removes 
completely the FL statate of YbInCu4. The experiments give [/, e ffh® = T v , where /i e // is the effective magnetization 
of the 4/ state in the LM regime, i.e., at the VF transition, the magnetic energy of the paramagnetic / states 
is comparable to the Kondo energy of the FL phase. The critical field is temperature dependent and the phase 
boundary between the LM phase and the FL phase is given by the expression h c (T) = h®yjl — (T/T v ) 2 , where /i° is 
the zero-temperature value. 

The aforementioned behavior of YbInCu4 and similar systems can be understood from the slave boson solution of 
the periodic Anderson model. Taking the chemical potential of the high-temperature phase close to the pseudo-gap 
of the unrenormalized c DOS we find a very smal Kondo temperature (see Sec. Ill C[) , which explains the Curie-like 
susceptibility, a large and nearly temperature-independent resistivity, small thermopower, and a negligible magnetore- 
sistance observed for T >T V . Since the Kondo screening is absent in the LM phase, we can neglect the hybridization 
altogether and discuss the high-temperature properties of YbInCu4 by adding the Falicov-Kimball term to the effective 
Hamiltonian. This term opens a gap or a pseudogap in the excitation spectrum 4 ^ and explains most of the qualitative 
features in a self-consistent way. For a quantitative description of the paramagnetic phase, one would have to include 
the corrections due to the CF splitting of the / states^. 

The peculiar feature of YbInCu4 and similar systems is that the local moments remain unscreened down to very 
low temperatures. System with a pseudogap in the c DOS cannot remove the paramagnetic entropy by Kondo effect 
but have to approach the ground state by an entirely different route. In YbInCu4 the transition into a low-etropy 
state is achieved by an iso-structural phase transition which expands the lattice and facilitates the valence fluctuations 
between xthe low-volume 4/ 13 and the large-volume 4/ 14 state of Yb. The presence of the 4/ 14 configuration in the 
ground state means that some / holes are transferred in the c band, so that the chemical potential is shifted away 
from the pseudogap. This increases the Kondo coupling and makes the Kondo temperature of the expanded lattice 
comparable to T„. For T < T%, the local moment disappears, because the / and c states form hybridized bands, 
i.e., the / holes participate in the Fermi volume. The f—c charge transfer induced by hybridization stabilizes the 
low-entropy FL state and compensates the loss of elastic energy due to the lattice expansion. Once the Kondo scale of 
the hybridized system is equal to T„, the expansion terminates, as can be seen from the following argument. The VF 
transition takes place at T„, where the free energy of the high-temperature phase, !Flm, is equal to the free energy 
of the low-temperature (expanded) phase, Tfl- Below T„, the free energy is dominated by the hybridization term 
in the internal energy which stabilizes the FL state, Tfl < J~lm- Above T„, where the free energy is dominated by 
the paramagnetic entropy, the local moments of the pseudogapped phase are completely free, while in the expanded 
lattice they would be partially screened. Since the free energy of the expanded lattice would be bigger than J~lm ■, the 
most stable configuration above T„ is the low-volume (paramagnetic) one. 

In the expanded lattice, the chemical potential is still rather close to the minimum of the c DOS, which makes the FL 
scale much larger than the Kondo scale T„ (see Sec. Ill D|) . If we estimate the FL corrections to the T — value of the 
magnetic susceptibility or the electric resistivity up to O[(T/T ) 2 ) and assume Tq/Tk < 10 (as indicated by the data), 
the maximum relative deviation at T ~ T v is 1%. The Kadowaki- Woods ratio in the FL phase is anomalously small, 
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because the CF splitting does not affect the delocalized / states, so that the / states are effectively 8-fold degenerate. 
The FL analysis presented in Sec. III. 3, shows that the A coefficient is then reduced by a factor 1/N(N — 1). Thus, 
the slave boson theory provides an overall description of YbInCu4-like compounds at ambient pressure. 

To account for the effects of pressure, we recall that pressure stabilizes the low- volume (4/ 13 ) configuration with 
respect to the large- volume (4/ 14 ) one and model this effect by shifting the / level away from the chemical potential. 
This reduces the f—c coupling and the Kondo scale, shifts T v to lower temperatures and removes eventually the VF 
transition. For very large pressures, the paramagnetic entropy of the LM phase is not removed by the the lattice 
expansion but by a transition into a magnetically ordered (MO) phase which takes place at Curie temperature T c . 

In the presence of an external magnetic field, the condition for the phase boundary, Tlm^T-,^) = Ffl(T, h), can 
be approximated by the condition SLM(T v ,h) — const., where Slm{T v1 K) is the entropy of the LM phased. This 
approximation uses the fact that T is not affected by the pseudo-gap, so that T 3> T v , and the temperature and the 
field-dependence of Tfl can be neglected for T <T v (h). The critical line obtained in such a way agrees very well with 
the experimental dat a 49 ' 51 and with the slave boson result given by Eq. ([39]). For large fields, the slave boson solution 
shown by the dashed-doted curve in Fig [5] explains the metamagnetic transition that takes place at the critical field, 
h° c ocT v . 

Similar reasoning explains also the unusual behavior of the newly reported system Yb3Pt 4 52 ' 53 . In the paramagnetic 
phase YbaPt4 has the Curie-like susceptibility and nearly constant resistivity which exceeds the spin-disorder limit. 
Like in YbInCu4, the LM phase of YbaPt4 looks semimetallic, with unrenormalized ('light') Bloch sates and localized 
/ states. In the absence of the f—c exchange coupling, the paramagnetic entropy cannot be reduced by the Kondo 
screening but the ground state has to be approached by an alternative route. In Yb 3 Pt 4 the ground state is magnetic 
and we speculate that the LRO is due to a spin density wave (SDW). The transition between the LM and the SDW 
phase takes place at T/v ~ 2 K. The data show that the effective mass of metallic electrons becomes heavy at T/v^ 3 -. 

This behavior can be explained by the slave boson approach, assuming that the chemical potential is close to the 
pseudo gap of the unperturbed c band and admitting an AFM solution. The generalized slave boson equations which 
include the coupling to the staggered magnetization provide at Tn a non-trivial solution with hybridized states and 
a large FS. The reason is that large internal fields shift the chemical potential of up and down spins away from 
the singularity in c DOS (like in the case of an external magnetic field, discussed in section HlF[) . Since the Kondo 
coupling in the LM phase is very small, the c electrons are nearly free and the Kondo screening can only occur at 
extremely low temperatures. Unlike YblnCiu, where the paramagnetic entropy is removed by an iso-structural phase 
transition which enables hybridization and the transition to the FL phase, in YbaPt4 the paramagnetic entropy is 
removed by a spin density wave (SDW) transition which partially gaps the Fermi surface. The formation of the 
SDW at temperature Tn ~2K switches-on the hybridization, which delocalizes the / states, changes the effective 
degeneracy of the / states, and renormalizes the metallic mass. The chemical potential in the SDW phase is still 
rather close to the pseudo-gap, so that the ensuing FL scale is large, T ^> T/v (see discussion in Sec. Ill D[) . Hence, 
the mass enhancement at the transition is moderate or small. It would be interesting to check the assumption about 
the pseudogap in the c DOS by performing the band structure calculations that would take into account the large 
Coulomb repulsion between the / electrons. 

3. The T ~ T K case 

As a final example we consider the anomalous pressure dependence of the T 2 coefficient of the resistivity, A(P), 
and the residual resistance, po(P), observed in heavy fermions like CeC^Si^, CeC^Ge? 5 - 5 - or CePd2Ge2 5 -£ at very 
low temperatures. In these compounds, the FL scale is about the same as the Kondo scale inferred from the low- 
temperature peak in p(T) or a(T) 17 i 55 . The scales To and Tk are much smaller than the CF splitting (estimated from 
the high-temperature peak in in p(T) or a(T)i£), so that the excited CF sates can be neglected for T <Tr- <C Acf- 
At ambient pressure, the ground state of these compounds is often superconducting or magnetic and the FL behavior 
is enforced by applying an initial pressure Pq. For P > Po, the data show that A(P) decreases gradually from large 
initial values, drops at a critical pressure P c by nearly two orders of magnitud o i 55 ' 56 , and then continues a gradual 
decrease. The residual resistance increases from a small initial value, rises rapidly to a sharp peak at P c , and decreases 
at very high pressure to rather small value s 54 ! 55 . 

We explain these features by assuming that for Po < P < P c the effective degeneracy of the model is defined by 
the lowest CF state and that pressure changes the f—c coupling and Tk but not Acf- The FL state established 
at Po has a large Fermi surface, because the / electrons are delocalized and participating in the Fermi sea. Taking, 
for simplicity, the half-filled conduction band and describing the 4/ state be a CF doublet, we find that the Fermi 
surface is close to the edge of the Brillouin zone and that the average Fermi velocity in Eq. (|4"2"|) is small. In the 
aforementioned compounds the FL scale is small, so that A(Po) large. As long as pressure does not change the 
degeneracy of the / states, Luttinger theorem preserves the Fermi surface, so that vp and A/u in Eq. (f?5] ) remain 
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approximately constant. Since TV is also constant in this pressure range, the pre-factor of the (7T) 2 -term in Eq. (|42|) 
does not change much. Thus, the main effect of pressure is a gradual reduction of A{P) from the maximum value 
attained at Pq- This reduction is due to the increase of Tq(P), as can be seen from the linear scaling between \J A{P) 
and the inverse Kondo scale of the system 1/T ^-(P)**. 

At large enough pressure the hybridization becomes too large for the system to support the CF excitations, so that 
the effective degeneracy of the ground state increases. For P > P c , the degeneracy is not set by the lowest CF level 
but by the full multiplet (a sextet, in the case of Cerium). The slave boson solution of the periodic Anderson model 
with the SU(N) symmetry and infinite correlation shows that the Fermi volume decreases as TV increases, because 
a single / electron has to be distributed over more and more channels. The chemical shift A/x and the specific heat 
coefficient 7 also decrease with TV, while v 2 F becomes larger on a smaller FS, so that A(P) drops sharply for P > P c . 
The degeneracy of the / state does not change for P > P c and an increase of pressure above P c reduces A(P) by 
increasing To. In this pressure range, we find the scaling between A(P) and 1/T K , where T 2 K is the Kondo scale of 
the TV-fold degenerate model. 

The residual resistance also tracks the pressure-induced changes in the effective degeneracy. If the / states are 
localized at ambient pressure, as seems to be the case with CeC^Ge^ or CePd2Ge2^ at very low temperatures, the 
initial values of po(Po) are small. Taking, for simplicity, the model with a ground state doublet and an excited CF 
quartet, and using the FL laws, we find that the temperature dependence of the resistivity at ambient pressure is due 
to two resonant channels, while most of the current is carried by four non-resonant channels which have temperature- 
independent conductivity. The contribution of the resonant channels to the residual resistance can be neglected. At 
large enough pressure, P ~ P c , the CF excitations are removed, the degeneracy of the ground state changes from 
doublet to sextet, and all the channels become resonant. Because of Luttinger theorem, the Fermi volume shrinks 
for two of the channels (former resonant ones) and expands for four of them (former non-resonant ones). Hence, for 
P ~ P c , the overall contribution to the residual resistance increases sharply. A further increase of pressure does not 
change the degeneracy of the model but gives rise to the charge transfer from the / to the c states. This reduces 
the residual resistance and gives the po(P) curve its asymmetric shape. In other words, the FL liquid laws and the 
slave boson solution of the periodic Anderson model show that the large peaks in A(P) and po(P) are due to the 
pressure-induced change in the degeneracy of the / state. 

IV. SUMMARY AND CONCLUSIONS 

It is well known that the degeneracies and the splittings of the 4/ states have a strong impact on the Kondo scale 
and the high-temperature behavior of intermetallic compounds with Ce, Eu or Yb ions. In this paper, we have shown 
that the details of the conduction electron band structure have a large impact on the ratio of the Kondo scale to the 
Fermi liquid scale, which determines the type of the crossover between the incoherent and coherent regimes. 

Our analysis is based on the periodic Anderson model, which can be mapped at high temperatures on a single 
impurity Anderson or Kondo models with Kondo scale Tk ■ The incoherent properties of the lattice can be understood 
using a single impurity approximation which takes into account the structure of the / states. At low temperatures, 
the slave boson solution of the periodic model yields the FL laws characterized by an energy scale To. The Kondo 
and the FL scales depend on the shape of the c DOS in the vicinity of the chemical potential, the degeneracy and 
the CF splitting of the / states, the number of c and / electrons, and their coupling. The ratio To/T^ is determined 
by the details of the band structure. Depending on the relative magnitude of T and Tk, the crossover between the 
high- and low-temperature regimes proceeds along very different routes. A sharp peak in the c DOS yields To ^CTk 
and gives rise to a 'slow crossover', as observed in YbAl3 and similar compounds. The minimum in the c DOS yields 
To ^>Tk, which causes an abrupt transition between the high- and low-temperature regimes, as observed in YbInCu4- 
like systems. In the case of CeCu2Ge2 and CeCu2Si2, where To —Tk, our results show that the pressure-dependence 
of the A(P) coefficient and the residual resistance can be related to the change in the degeneracy of the / states. 

The slave boson solution shows that the low-temperature response of the periodic model can be enhanced (or re- 
duced) with respect to the predictions of the single- impurity model with the same Kondo scale. In the coherent regime, 
the renormalization of transport coefficients modifies the Wiedemann-Franz law and can lead to an enhancement of 
the thermoelectric figure-of- merit. The FL laws explain the correlation between the specific heat coefficient 7 and the 
slope of the thermopower, a(T)/T, or between 7 and the T 2 coefficient of the electrical resistance A = p(T)/T 2 . In 
the case of a TV-fold degenerate model, the FL laws explain the deviations of the Kadowaki- Woods ratio Rrw = A/j 2 
and the q ratio, q = \e\lim^x^o}x a /"fT , from the universal values. 

The shape of the c DOS affects also the magnetic response of the system. The field-dependence of the magnetization 
resembles the temperature-dependence of the entropy, as can be seen from Figs.[5]and O The c DOS with a maximum 
around the chemical potential leads to an extended plateau in m(h), while a pseudo-gap in the c DOS leads to a 
metamagnetic transition for the fields of the order of the Kondo temperature. 
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APPENDIX A: GENERALIZATION TO MULTI-ORBITAL SYSTEMS 



1. General formalism 



Here, we generalize our model to more realistic systems where the local f orbital have a supplementary degeneracy. 
The later is lift at low temperature by the crystal field. In the limit U — > oo, the PAM Hamiltonian (JTJ is generalized 
to 



ktr 



(Al) 



where a — 1, • • • ,N/2 is a local orbital index, and Ag F is the crystal field. For Ce and Yb we have N = 6 and 
N = 8. The slave boson approach described in section [TT] for N = 2 can be easily generalized within the mapping 
\0)( — ► b\ |0)j, \aa)( — ► d\ aa |0}., with the exclusion of all the states with more than one electron on a given site. 
This leads to the local identities b\bi + J2aa fiavfiaa — 1) which are satisfied by introducing Lagrange multipliers Aj. 
Within the mean-field approximation, we replace the bosonic fields bi and the Lagrange multipliers Aj by homogeneous 
static real fields r, and A. The quadratic mean-field Hamiltonian © obtained for N = 2 is then generalized to 



Ef-X-fj, 



(l-r 2 ) 



(A2) 



The self-consistent parameters r and A are obtained by minimizing the free energy (3J-(r, A) = — In Tre @ H , such that 
dT{r, X)/dr = and dT(r, X)/dX = 0. This gives 



2r(E f -X-fi) = ^^(cLd kQCT + 4 QCT c kff } , 

kacr 
kao" 

where, (• • •) denotes the thermal average with respect to the mean-field Hamiltonian (|A2[) and 



(A3) 
(A4) 

(A5) 



The average electronic occupation per site of the / and c orbitals is nj = ^2 i(J {f} ai yfiaa) — J2ia(dia<jdiaa) and 

n c = X)io-( c L c io-)' respectively. These averages are determined by the slave-boson amplitude and the total electronic 
occupation, 



5>7 - 1 



n r = n 



The self-consistent Eqs. (|A3[ IA4I IA5[1 can be rewriten as 

Ef-X-n 



J —oo 

Pd(u)n F (uj)duj , 



T 



+00 



p c {io)nF{oj)duj 



(A6) 
(A7) 

(A8) 
(A9) 
(A10) 
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where np(ui) = l/[l + e^ w ] is the Fermi function and pdc, Pd an d p c are the spectral densities of the local single-particle 
Green's functions of a given spin. These Green's functions are defined in the usual way as thermal averages of the 
(imaginary) time-ordered products of the appropriate creation and annihilation operators. 



Pd 



H = --£E ImG ^( k > 

7T * — ' » — ' 

k a 

k 

(w) = --Vim G c (k,cj) . 

7T fc * 



p. 



The Green functions are computed from the quadratic Hamiltonian (|A2|) . We find 

1 



G cc (k, w) 



w + H - £k - r 2 ^ 2 £ Q V(w - A - Ag F ) ' 
rV 



lu - X - A 



— G cc (k,w) 



u - X - Ag F (w - A - Ag F ) 2 



G cc (k, w) 



We define the non-interacting electron Green's function 

GoH=x;— 



W - £k 



The local spectral densities can be rewritten as 

Pc{uj) 



7T ^ „ a; — A — Ag F y 



i. 



-Im 



-— Im 

7T 



A - A 



w — A — Ag F (a; — A — A CF 



w - A - Ag F 



(All) 
(A12) 
(A13) 

(A14) 
(A15) 
(A16) 

(A17) 

(A18) 
(A19) 
(A20) 



2. Degenerate case: no crystal field splitting 

We consider here A^ F = 0. In this case, we can check easily that the complete formalism described in this article 
for N — 2 can be generalized to any value of N, within the following rescaling 



V 2 -> *V 2 



(A21) 



Note that, due to the exponential dependence of the energy scales To and Tk with respect to V 2 (see Eqs. (|20j|22ll54"l) . 
these two scales have to be rescaled like 



T {N) _ T (2) -(2-A0/ATJ( Mo )po(po) 
1 K — 1 K e 



T (N) _ T (2) -(2-N)/NJ( Jl0 )p (ji ) 

where (N) refers to the number of degenerate channels considered for the / orbitals. 



(A22) 
(A23) 



3. Kondo temperature 

For r — ► 0, the generalized mean-field equations can be written as 

n c f + °° 

y = / po{uj + p,)n F (uj)duj , 



(A24) 
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2 



= n F (\ + A% F ) , (A25) 
= E ,- A -^ + M) tanh (_) ^ . (A26) 



We can arbitrarily choose one vanishing Ag,^, = 0, corresponding to the doublet with the lowest local energy Ef. 
If, for the other doublets, one has >> T, then we obtain n™^° = 0, and the formalism developed for N = 2 

can be applied simply. In the more general case where A^° ~ T, or A^° < T then one has to consider also the 
contribution from these orbitals. 
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